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■ IT IT IT correlations from Pb+Pb collisions at 158 GeV/c per nucleon are pre- 
sented as measured by the focusing spectrometer of the NA44 experiment at CERN. 
The three-body effect is found to be stronger for PbPb than for SPb. The two- 
dimensional three-particle correlation function is also measured and the longitudinal 
extension of the source is larger than the transverse extension. 
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1 INTRODUCTION 

Three-particle Bose-Einstein correlations are sensitive to source geometry and chao- 
ticity[^ in a way not seen in ordinary two-particle correlations. Hence three-body 
correlation studies provide a test of the validity of various assumptions often used in the 
parameterization of the two-particle correlation function. The additional information is 
carried by the genuine three-particle correlation term. In the case of source asymmetries 
and/or coherent particle emission, the strength of the genuine three-particle correlation 
term will be suppressed. Genuine three-particle correlations have been found at their full 
strength in electron-positron annihilations and in high energy in pp collisions[0. In 
other reactions and at lower energies, the possible presence of a three-body correlation is 
unclear 1^, |], |^, ||. In heavy-ion collisions a significant suppression of the genuine three- 
body correlation has been reported by this collaboration in S+Pb collisions^. In this 
paper we have analysed our Pb+Pb three-pion events in the same way as in . In addition 
we report on the first experimental measurement of the two-dimensional three-particle 
correlation function. 



2 EXPERIMENTAL SET-UP 

The NA44 experiment is a focusing spectrometer measuring particle distributions 
at mid-rapidity with excellent particle identification. The spectrometer set-up has been 
described well in [^. The momentum range and sign of the charge of the particles is 
selected by two dipole magnets. The data used for this analysis is 4 GeV/c ± 20%. The 
spectrometer axis is located at 44 mrad with respect to the beam axis and covers a pt 
range of - 400 MeV/c, with an average of (pr) = 145 MeV/c. The rapidity range is 
3.1 — 4.1, with an average of 3.7. The trigger requires a well identified single lead ion, at 
least two hits on the hodoscopes and a pion signal (with no electrons) from the Cherenkov 
detectors. The centrality is fixed at the most central 9% of the geometrical cross-section 
by means of a threshold on a scintillator downstream of the target. 



3 THE PARAMETERIZATION OF THE CORRELATION FUNCTIONS 

In the case of a totally chaotic source, the two- and three-particle Bose-Einstein 
correlation functions can be written as0: 

C2 = l + |F,,f, (1) 
Cs = 1 + |Fi2|' + |F23|' + |F3i|2 + 2i?e(Fi2F23F3i), (2) 

where, assuming plane wave propagation: 

F,, = J e^«->(r)rfV, tj = 12, 23, 31 (3) 

where Qij is the four-momentum difference of particle i and j and p is the source density 
function. The last term in Eq. (Q) is the so-called genuine three-particle correlation due 
to a pure three-body effect. By assuming a symmetric Gaussian source density function 
p of width R the correlation functions are parametrized as: 

C2{Q^,) = l + Ae-«-^^', (4) 
Cs{Qu,Q23,Q3i) = l + AEe-^?'^'+2A3/V^(^•^«?.)^^ (5) 
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where ij is one of the permutations ij = 12,23,31. Here A is a phenomenological ||1 1|| 
parameter defined by: A = C(Q = 0) — 1. This parameter has been introduced due to the 
fact that the measured two- and three-particle correlation functions do not reach the full 
value which, for Bose-particles, are 2 and 6 for two and three particles respectively. 

We have analysed the three-particle correlation function in terms of the sum of the 
momentum differences, Q3, and the components of Q3 transverse and along the beam, i.e. 
Ql = Qt + Qf- We will use Eqs. (^, (||), and (|^) to parameterize our data: 

C^3(g3) = l + A3e-«i^i, (6) 

where Q| = QI2 + Q23 + Qii^ and: 

C^3(g*,gO = l + A2.e-(«?Q*+^?«?), (7) 

calculated in the longitudinal center-of-mass system for the triplets characterized by 
J2Pz = 0- The kinematical variables used in Eq. (J^) are defined in Eqs. 



Qt — Y Q't,12 + '?T,23 + 5't,31; 'iT,ij — iPx,i Px,jY + {Py,i PyjY 



Ql = y(lL,12 + 9l,23 + ?i,31. (ihj = iPz,i - Pzjf (9) 

where px-, Py, and pz are the three momenta components. 

4 THE GENUINE THREE-PARTICLE CORRELATION 

We introduce a weight factor uj defined by the relation: 

C3 = 1 + + + |F3i|2 + 2IF12IIF23IIF31I X LU. (10) 

The weight factor a; is a measure of the strength of the genuine three-particle correlation. 
It can be experimentally determined by using the following expression, extracted from 
Eq. (0): 

{CsiQs) - 1} - {C2{Ql2) - 1} - {C2{Q23) - 1} - {C2{Q3l) - 1} 



2^J{C2{Ql2) - 1}{C^2(Q23) - l}{C2{Q3l) " 1} 

For a totally chaotic and symmetric source u = 1, but u will differ from 1 for an asym- 
metric and/or coherent particle emitting source, see 0, ^ H, |12], |l^. 



Traditionally, the chaoticity is measured by the intercept parameter A of the two- 
particle correlation function. In practice this intercept measurement is difficult since phase 
space approaches as Qij — > and C2 is diluted by resonances and smeared out by the 
effects of momentum resolution and coulomb repulsion. The u is, however, calculated over 
a relatively broad range in Q. Any deviation of u from 1 is an indication of coherence 
as asymmetry effects have only minor influence on Note that lambda cancels in Eq. 
(|TT|) . This is seen directly for Qij = as a;(Q = 0) = 1 for all values of lambda. 



5 DATA ANALYSIS 

We study negative pions produced when a lead beam hits a lead target. The dataset 
consists of 94,000 events with three pions. We flnd tracks by fltting straight lines to hits on 
two hodoscopes, two strip chambers and a pad chamber, all situated behind the magnets. 
Tracks are not allowed to share the same hodoscope slats and a minimum track separation 
is required in the pad-chamber to certify a high purity of the three-pion sample. The 



time-of-flight start signal is derived from a beam counter with a time resolution of a ^ 35 



ps|jTj]. Events containing particles other than pions are rejected by combining information 



from the time-of-flight and from the multi-particle threshold imaging Cherenkov (TIC) 



detector ll^]. The TIC distinguishes pions from heavier particles on a track by track basis 
by requiring a sufficiently high signal in a fiducial zone around the tracks. The residual 
contamination from particles other than pions is typically less than 1%. Identical cuts are 
used in constructing the two- and three-particle correlation functions. 

The bin size used in the analysis corresponds to the momentum resolution. The 
resolutions are calculated using a Monte-Carlo programme with full detector simulation. 

The correlation function is determined using C{Q) = A{Q)/B{Q). The "real" mo- 
mentum distribution A{Q) is constructed from tracks from the same event and the "back- 
ground" distribution, B{Q), is constructed from tracks mixed randomly from all events 
contained in A{Q). Around ten "background" events are created for each "real" event 
in order to avoid statistical uncertainties from the "background" sample. This method 
cancels effects of the experimental acceptance and trigger biases and is described in detail 
in our previous publication |T^. 

To compare with theoretical correlation functions, corrections are applied iteratively 
to produce the correlation function Ccorr- 

Ccorr Craw ^ Kgp(j X Kacceptance ^ K(joulomb' (-^^) 

This procedure converges within four iterations. The factors are explained below: 

• The background spectrum is distorted with respect to the true uncorrelated many- 
particle spectrum, owing to the effect of the many-particle correlations on the single- 
particle spectrum. This is iteratively corrected by the factor Kspc appropriately gener- 
alized to the three-particle case, in which each particle used in the background spectrum 
is weighted by the correlation from the event from which it is taken. 

• The factor Kacceptance corrccts the data for the momentum resolution of the spec- 
trometer and the many-particle acceptance and is calculated using a Monte Carlo pro- 
gramme with a full simulation of the tracking detectors and multiple scattering. 

• The two-particle correlation function has been corrected for the Coulomb in- 
teraction, Kcouiomb, between the particles using the Coulomb wave function integra- 
tion method [p!7|, IT8[. The three-particle correlation function has been corrected for the 



Coulomb interaction using a similar technique |jT9|, ^ by using a three-body Coulomb 
wave-function pO, 511. 



Coulomb interactions with the residual nuclear system are expected to be small 
and no corrections are applied in this analysis. Final-state strong interactions are also 
expected to be small and due to large uncertainties in proposed procedures, no corrections 



are applied for them 122 



The systematic errors are evaluated by varying the analysis parameters and cal- 
culating the difference in the correlation functions produced. These variations include 
changing the momentum resolution assumed in the Monte-Carlo correction by ±10%, 
changing the time-of-flight cuts, increasing the minimum separation in the pad chamber, 
requiring more hits in the strip chambers, and increasing the minimum slat separation 
in all the hodoscopes. We re-calculate u for each new set of correlation functions. The 
systematic errors are estimated by summing up the differences to the mean-value for each 
altered setting. The statistical error on the u distribution is: astati'^) = ct/VN, where 
a is the variance and N is the number of entries. The statistical error is much smaller 
than the systematic error. In the statistical error calculation we have not included error 
propagation from the correlation functions to uj. This is very complicated so instead we 



As 


Rs (fm) 




1.92 ± 0.16 ± 0.33 


3.78 ± 0.16 ± 0.37 


17.1/15 



Table 1: Results from a fit to the three-pion correlation function using Eq. (|^). The errors 
are statistical and systematic respectively. 



System 


A 


R (fm) 




Centrality 


37r ^ 27r 


0.57 ± 0.04 


7.49 ± 0.34 


18/20 


9% 


2vri 


0.52 ± 0.04 


7.56 ± 0.38 


30/36 


18% 



Table 2: Results from a fit to the two-pion correlation function made from the three- and 
a two-particle dataset. Both datasets are corrected by the Coulomb- Wave integration but 
they have different centrality. The errors are statistical only. 



have calculated u using Cs±crstat and C2±crstat and treated the deviation as an additional 
systematic error. The two systematic errors are added in quadrature. 

There is still some room for further systematic errors in the Coulomb correction, 
due to the fact that the unknown exact three-body Coulomb wave-function is reproduced 
only asymptotically. However, in the kinematic region of the NA44 experiment, these 
non-asymptotic correction terms in the three-body Coulomb wave-function are known to 
decrease strongly with increasing energy of the triplets and the average energy of the 
triplet is large compared to the typical scale of the three-body Coulomb potential in the 
NA44 three-pion data sample, see |jl9| . 



6 RESULTS AND DISCUSSION 

The three-pion correlation function is shown in Fig. |1] and the result of a fit to Eq. 
d^) is summarized in Table |]. An estimate of the strength of the three-pion correlation 
function is shown in Fig. ^ as dotted lines. This is done using Eq. (|^) with A and R 
extracted from a fit to the two-pion correlation function. The lower dotted line is without 
contribution from the genuine three-body correlation, i.e. the datapoints would follow this 
path in case of uj=0. The upper dotted line is for u =1, and in fact seems to follow the 
datapoints. 

The two-pion correlation function is constructed by using the three combinations of 
a pair from the triplet datasample. The correlation function is fitted using Eq. (Q) and the 
results are summarized in Table ^ where we also list our results from . The extracted 
fit-parameters are consistent between the two- and three-particle datasets despite the 
different centrality. The similarity indicates that the effect on the correlation function due 
to the presence of a third pion is small for our data. By comparing results in Table ^ and 
ID one sees directly that A3 > 3A as expected when u > 0. 

The direct method of extracting the strength of the genuine three-pion correlation is 
to calculate the u factor using Eq. (0). When we have determined C2 and C3, as described 
above, we can use the actual data points, C2{Qi2), C2{Q2z)i C'2(Q3i), and C^iQ^), so as 
not to be biased by some parameterization. The data points are obtained by using Q12, 
Q23) Q31, and Qs for each event and we have checked that the result do not change when 
altering the Q's by ± 5%. In order to avoid poles in the denominator of Eq. ([TT|) events 
are accepted if Qij < 60 MeV/c. As a result we obtain a distribution of u for each 
bin, see Fig. ^, where we also show our S-I-Pb result 0. Events are taken in the region of 
the genuine three-body correlation, i.e. in the range 14 < Q3 < 54 MeV/c as we have no 
data below 14 MeV/c. In this range we find the weighted-mean uJ = 0.85 ± 0.02 ± 0.21 
where the errors are statistical and systematic, respectively. The value uJ from our S+Pb 



A 



^2d 


Rt (fm) 


Ri (fm) 




2.67 ± 0.49 ± 0.19 


3.67 ± 0.26 ± 0.19 


5.89 ± 0.53 ± 0.68 


189/113 



Table 3: Results of a fit to the two-dimensional three-pion correlation function using Eq. 
(0). The errors are statistical and systematic respectively. 



result was uJ = 0.20 ± 0.02 ± 0.19. 

The three-pion correlation in PbPb collisions has also been published by WA98p3|| 
without using the more correct Coulomb correction method[^ employed by NA44. Within 
the errors, however, the uJ are compatible between the two experiments. 

The two-dimensional three-pion correlation function has also been analyzed in the 
longitudinal-center-of-mass system. Only bins which contain more than 40 entries are used 
in the fit. Projections are made following and are shown in Fig. ^ The results of a 
fit to the two-dimensional correlation function using Eq. (J^) is summarized in Table |^. 
The systematic error of Rt and Ri are correlated, i.e. the radii both become small or both 
large when varying the analysis parameters. When comparing the difference of the radii 
to the largest statistical error we see a 4-a effect in the difference of the radii. 

There exist only few predictions of the two-dimensional three-body correlation func- 
tion. One is based on the Lund string- model [p5[] where a larger longitudinal than transverse 
radius parameter is in fact expected. 



7 CONCLUSIONS 

We have found that the measure of the strength of the genuine three-particle cor- 
relation, expressed as the mean weight factor u in PbPb interactions, is compatible with 
1, 00= 0.85 ± 0.02 ± 0.21. This is different from SPb interactions where we earlier have 
found, using the same analysis method, C(j=0.20 ± 0.02 ± 0.19. We consider this differ- 
ence to be significant as the systematic errors between the two dataset are correlated. 
The small a;-value in SPb interactions indicates a nonchaotic mechanism for particle 
production[^ ^ ^ different from PbPb interactions which are compatible with a 
fully chaotic mechanism. 

We have performed a first measurement of the two-dimensional three-pion correla- 
tion function. The data show that the longitudinal radius parameter is larger than the 
transversal radius parameter in the longitudinal center-of-mass system. 
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Figure 1: The three-particle correlation 
function fitted (solid line) using Eq. (^. 
The double error bars are statistical (in- 
ner) and systematic (outer) respectively. 
The dotted lines indicate what path the 
datapoints should follow with (upper) and 
without (lower) the genuine three-particle 
correlation, see the text. 



3.5 









I 












■ '1 1 


□ Pb+Pb 
• S+Pb 




1 1 1 



20 40 60 



Q3(MeV/c) 

Figure 2: The weight factor a; as a func- 
tion of Qs for Pb+Pb (open box) and 
S+Pb (closed circles). The dashed area in- 
dicates the Qs-range without data. The er- 
rors shown are systematic only since the 
statistical ones are within the data points. 
Bear in mind that the sytematic errors are 
correlated within each data set as well as 
for the two data sets (see text). 
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Figure 3: The two-dimensional three-pion correlation function is projected onto the A) 
Qi and B) Qt axis for a range in Qt and Qi as indicated in the figure. The solid line is 
drawn using Eq. (0) using an average {Q) in the orthorgonal direction of the momentum 
difference. The correlation function is calculated in the longitudinal center-of-mass system 
for the triplets. The errors are statistical (inner) and systematic (outer) respectively. 
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